INVERSE OBSTACLE PROBLEM FOR THE NON-STATIONARY 
WAVE EQUATION WITH AN UNKNOWN BACKGROUND 
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Abstract. We consider boundary measurements for the wave equation on a 
bounded domain M C M 2 or on a compact Riemannian surface, and intro- 
duce a method to locate a discontinuity in the wave speed. Assuming that the 
wave speed consist of an inclusion in a known smooth background, the method 
can determine the distance from any boundary point to the inclusion. In the 
case of a known constant background wave speed, the method reconstructs a 
set contained in the convex hull of the inclusion and containing the inclusion. 
Even if the background wave speed is unknown, the method can reconstruct 
the distance from each boundary point to the inclusion assuming that the Rie- 
mannian metric tensor determined by the wave speed gives simple geometry 
in M. The method is based on reconstruction of volumes of domains of influ- 
ence by solving a sequence of linear equations. For r £ C(dM) the domain 
of influence M(t) is the set of those points on the manifold from which the 
distance to some boundary point x is less than t(x). 



1. Introduction and the statement of the results 
Let us consider the wave equation on a compact set Mcl 2 , 
(1) d 2 u{t, x) - c(x) 2 Au(t, x) = 0, (t, x) £ (0, oo) x M, 

u(0, x) = 0, d t u(0, x)=0, xe M, 

d v u(t, x) = f(t, x), (t, x) e (0, oo) x DM, 

with a piecewise smooth wave speed 




cq{x) + cs(x), 

cq(x), x € M \ £, 



where £ C M mt and Co and cs are smooth strictly positive functions. Here d v 
denotes the normal derivative and the boundaries dM and 9E are smooth. Physi- 
cally £ corresponds to an obstacle or an inclusion in which acoustic waves propagate 
faster than in the background medium modelled by cq- 

Let E and cs be unknown and let us assume either that the wave speed cq is 
known or that it is unknown and gives simpl^] geometry in M, We describe a 
method to locate the inclusion E using the operator 

(2) A 2T :/.W| (0 , 2T)x eM, feL 2 ((0,2T)xdM), 
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where u^(t, x) = u(t, x) is the solution of and T > is large enough. The oper- 
ator A 2 t models boundary measurements and is called the Neumann-to-Dirichlet 
operator. The method works equally well with an anisotropic background wave 
speed, see the wave equation ^ below. 

In a recent article [13] . Chen, Haddar, Lechleiter and Monk introduce a method 
to reconstruct an impenetrable obstacle in the Euclidean background by using 
time domain acoustic measuments. Their method is similar to the linear sam- 
pling method originally developed for an inverse obstacle scattering problem in 
the frequency domain |14j . The method we propose here is based solely on control 
theoretic approach in the time domain. Our method is similar to the iterative time- 
reversal control method by Bingham, Kurylev, Lassas and Siltanen [7j [15] originally 
developed to reconstruct a smooth wave speed as a function. Computationally our 
method consists of solving a sequence of Tikhonov regularized linear equations on 
L 2 ((0, 2T) x dM), and allows for a very efficient implementation if computation 
steps are intertwined with measurement steps. 

By using the boundary control method, a smooth wave speed can be fully recon- 
structed from the Neumann-to-Dirichlet operator. This is a result by Belishev [3] 
for an isotropic wave speed and by Belishev and Kurylev [5] for an anisotropic wave 
speed. However, the boundary control method in its original form is exponentially 
unstable and hard to regularize. To our knowledge, there are two computational 
implementations of the method [4j[26]. In addition to these two implemetations, 
the only numerical results related to the boundary control method we are aware 
of are in the recent article by Pestov, Bolgova and Kazarina [39] . We believe that 
there is a demand for methods that reconstruct less but are more robust. 

Our method locates the inclusion E by computing the travel time distance from 
each boundary point to E. Let us next describe the travel time distance function in 
detail. We denote by (M, g) a smooth complete and connected Riemannian surface 
that models the background wave speed. For example, in the isotropic case we 
have M = M 2 and 

(dx 1 ) 2 + (dx 2 ) 2 - 

g(x):= K - ' { ' , xgR 2 . 

c (x) z 

We let M C M and E C M mt be compact sets with smooth boundary and non- 
empty interior. Moreover, we assume that M is connected and that the wave speed 
with the inclusion is given by the non-smooth Riemannian metric tensor, 

(3) = ( c( f '* kix) > X H XT *,* = !, 2, 

yg jk (x), xGM\T,, 

where c is a smooth scalar function on E satisfying c{x) > 1 for all x G E. We 
denote by d(x,y), x,y G M, the Riemannian distance function of (M, g) and by d 
and d the Riemannian distance functions of (M, g) and of (M,g), respectively. Note 
that d is not necessarily a restriction of d. For example, if g is the Euclidean metric 
tensor and M is non-convex then there exist x, y G dM such that d{x, y) < d{x, y). 

In the case of a known background manifold (M, g) without boundary, our 
method reconstructs the boundary distance hull of E, 

H 9M (X):=M\ |J B(y,d(y,E)), 

yedM 
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FIGURE 1. The domain M is the large disc and the background 
wave speed is constant. The inclusion E consist of the two small 
discs. On left. If the backgroud is known, the method can recon- 
struct the boundary distance hull HqmC^) visualized as the shaded 
area in the picture. On right. If the background is unknown and 
simple, the method can reconstruct the vectors fs (y) grad rs (y) , 
y G dM, visualized as the line segments in the picture. 

where B(y,r) := {x G M; d(x,y) < r} for y G M and r > 0. If (M,g) is simple, 
then the embedding M C M plays no role in our method, and we can replace d 
with d in the definition of Hqm{^)- Note that E C Hqm{^) and, in the special 
case of the Euclidean background, the boundary distance hull Hqm{^) is a subset 
of the convex hull of E. 

In the case of an unknown simple background manifold (M, g) with known 
boundary (dM, g\dAi), our method reconstructs the distance function 

r E (y):=%,S), y G dM. 

Thus it determines the set 

B aM (E) := {(y,n) G dTM; \n\ g < d(y,E)}, 
where dTM := {(y,rf) G TM; y G dM, r) G T y M). The set B dM (Z) satisfies 
(4) iWE) = M \ exp(S aM (E)), 

where exp is the exponential map of (M , g) . As g is unknown, the equation Q 
does not give a method to determine Hqm{^)- However, if we know a priori that 
g is close to a known metric tensor g° on M, we get a distorted image of Hqm(E) 
simply by visualizing the set exp £?,9m(E), where exp is the exponential map of 
(M,g°). _ 

Moreover, if M is also simple and <?| is known, then the method can recon- 
struct the distance function, 

?t,{x) := d(x, E), x G M \ M. 

From this function it is easy to extract also directional information, since grad fs (x) 
exists for almost all x G M \ M and is a unit vector pointing to such a point z G 9E 
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that 

d(x, z) = d(x, E). 

See Figure [I] for a visualization of Hqm(£) and the directions giad r^(y), y G dM. 
In this paper, we focus on the recovery of the distance function fjj and do not study 
further the directional information contained in fjj. 

It is well-known that the high frequency behavior of the scattering pattern de- 
termines a convex obstacle in the Euclidean background. For a review of this and 
related results we refer to the survey article [40] . We emphasize, however, that our 
method does not rely on analysis of high-frequency solutions. Moreover, it seems 
possible that a logarithmic type stability result for the proposed method could be 
proven by using a stability estimate for the hyperbolic unique continuation prin- 
ciple. Such an estimate is formulated in an unfinished manuscript by Tataru [461 
Thm. 3.45]. Furthermore, our method is similar to the iterative time-reversal 
control method, and that method can be modified to work in the presence of mea- 
surement noise [7J. Such robustness against noise is typically not possible for a 
high-frequency solutions based method. 

From a point of view of numerical computations, our method consists of recon- 
struction of volumes of domains of influence. We define for a function r : dM — > 
[0, oo ) the domain of influence with and without the inclusion, 

M(t) := {x € M; there is y € dM such that d(x,y) < r(y)}, 

M(t) := {x e M; there is y € dM such that d(x,y) < r(y)}. 

Moreover, we define 

C T {dM) := {t e C(dM); < r{x) < T for all x e dM}, 

and denote by m and m the Riemannian volume measures of (M,7f) and (M,g), 
respectively. Using the method introduced in |37j we can compute the volumes, 

(5) m(M(T)), t e C T (dM), 

from the operator A2T by solving a sequence of linear equations on the space 
i 2 ((0, 2T) x dM). We outline this method in section [H] where we also general- 
ize it to cover a wave speed given by a non-smooth Riemannian metric tensor of 
the form ((3). 

Note that, on the one hand M(r) grows as c{x) in {[3J) grows, but on the other 
hand the factor c(x)~ 2 in the volume measure m(x) = c(x)~~ 2 m(x) gets smaller. 
We show that the latter effect is dominating near the boundary <9E. That is, we 
show the following result in section [2] 

Theorem 1. Let t : dM — > E be such that 

(6) M(r) n E ^ 0, M(t) n E mt = 0. 

Suppose that there is a neighborhood Us of E such that M{t) n Us is an embedded 
smooth manifold with boundary. Then 

(7) hminf rn(M(r + e))- MM(r + e)) > Q 

If the background manifold (M, g) is known, y £ dM and the sets, 
B(y , r) := {x € M mt ; d(x, y ) < r}, r> 0, 
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have smooth boundaries, then Theorem [T] gives a test to determine the smallest 
r > such that B(y ,r) HE^I, that is, the distance d(yo,T,). Indeed, we may 
choose r r £ C(dM) such that r r (yo) = r and t t — outside a small neighborhood 
of yo in dM. If T r (y) decreases to zero fast enough as d(y, yo) grows, then M(r r ) = 
B(yo,r). Moreover, M(r r ) = M(r r ) whenever M(r r ) flE = 8, whence by Theorem 

E 

d(y , S) = max{r > 0; m(M(r r )) = m(M(r r ))}. 
In section [3] we present a refinement of this test that works for compact domains 
M C M with smooth boundary, where M is a complete smooth manifold without 
boundary. In particular, (M,g) may have conjugate points and B(yo,r) may have 
non-smooth boundary for some yo £ dM and r > 0. Moreover, we show that in the 
case of an unknown and simple background manifold (M,g), the distance d(yo,Ti) 
can be reconstructed by using a test related to the second derivative of the function 
r i-> fh(B(y ,r)). 

^.Let us next summarize our results as a theorem. We remind the reader^ that 
(M, g) is a smooth complete and connected Riemannian surface, M C M and 
S C M mt are compact sets with smooth boundaries and non-empty interiors, M 
is connected and g is defined by Q . We consider the wave equation 

(8) dfu(t,x) - Api(t,x) = 0, (t, x) £ (0, oo) x M, 

u\ t =o = 0, d t u\ t=0 = 0, 

d v u(t, x) = f(t, x), (t, x) £ (0, oo) x dM, 

where Ag is the Laplace-Beltrami operator on (M, ~g) and <9„ is the normal derivative 
on (dM,g). We denote by \g(x)\ the determinant of g(x). If (g^ k (x))^ k=1 denotes 
the inverse of g(x) = (gjk(x))j k=i m some coordinates, then 

A~ g u:= 1^)1"^ (WI^W^ 
j,k=i v 

d u u:= g jk {x)v k {x)^T, 

where v = [y\, v%) is the exterior co-normal vector of dM normalized with respect 
to g, that is, Yl\k=i9 : ' kv j v k = 1- We define the operator A2T by |2| where is 
the solution of (181. 



Definition 1. A compact Riemannian manifold (M,g) with boundary is simple 
if it is simply connected, any geodesic has no conjugate points and dM is strictly 
convex with respect to the metric g. 

Theorem 2. If the boundary (dM,g\gM) o,nd the operator A2T are known, then 
the volume data i5L) can be computed. Moreover, if T > \\Ts\\ L00 rg M y then the 
following two implications hold. 

(i) If(M,g) is simple, then the set BqmC^) C dTM can be reconstructed from 
the volume data |5|). Furthermore, the boundary distance hull satisfies 

H dM m =M\exp(B aM (S)). 

(ii) If (M, g) has no boundary and g is known, then the boundary distance hull 
Hqm(^) can be reconstructed from the volume data pi). 
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In the context of inverse acoustic scattering problems, there is an extensive lit- 
erature about inclusion detection methods. Some of these methods have also been 
applied to solve inverse obstacle problems for the time domain wave equation. The 
methods in [33] and in [3 6) take Fourier transforms of time domain measurement 
data and solve the inverse obstacle problem by using inclusion detection meth- 
ods developed for scattering problems. Moreover, the method in [12] process the 
measurement data partly in the frequency domain. 

The only inclusion detection method processing the measurement data entirely 
in the time domain that we are aware of is the already mentioned sampling method 
in |13) . The analysis in [13] depends on frequency domain techniques, and the finite 
speed of propagation for the wave equation seems to be an obstruction in carrying 
out the analysis. On the contrary, our method is based on the finite speed of 
propagation and the complementary unique continuation principle by Tataru [47] . 

Well-known inclusion detection methods in the frequency domain include the 
already mentioned linear sampling method by Colton and Kirsch and the enclosure 
method by Ikehata [23 . A modification of the linear sampling method by Kirsch 
is called the factorization method |29j , and it can be interpreted by using localized 
potentials pjj] . The enclosure method is the first inclusion detection method based 
on the complex geometrical optics solutions developed by Sylvester and Uhlmann in 
their fundamental paper [T5]. For later complex geometrical optics solutions based 
methods see [24l I2T1 l4"9j . 

The factorization method has been applied also to electrostatic measurements 
[T5I ITT) , and the enclosure method was developed for both acoustic scattering and 
electrostatic measurements from the very beginning. For other methods to solve 
inverse obstacle problems related to scattering and electrostatic measurements see 
the probe [22] and singular sources [UJ methods, the no response test [35], the 
scattering support techniques [3TJ|43l[T9] and the review article [42]. Furthermore, 
we refer to the review article [25] for uniqueness and stability results related to 
inverse obstacle problems. 

The uniqueness results for the inverse problem for the wave equation mentioned 
above assume smooth wave speed [3[ [5J. However, in a recent article [25], Kir- 
pichnikova and Kurylcv consider piecewise smooth wave speeds on Riemannian 
polyhedra. Moreover, the stability results [TJ [44] establish uniqueness for wave 
speeds with a limited number of derivatives, and there is an extensive literature 
about uniqueness results for the related Calderon's inverse problem under low reg- 
ularity assumptions including [2] [HI EH [T7], [301 [35] . For a review of the latter results 
we refer to 05] . 

2. Inclusion detection from the volume data 

In this section we prove Theorem [TJ We assume throughout the section that 
r : dM ->■ R satisfies @. 

Lemma 1. Let e > 0. If x E M(r + e) and 

(9) leS or x<£M(r + e), 

then there is z E such that 



(10) 



d(x, z) + d(z, M(t)) < e. 
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Proof. Let 1(7) denote the length of a path 7 with respect to metric d and £(7) with 
respect to metric d. There is y € dM and a path 7 : [0,£] — > M from y to x such 
that 7(7) < r(y) + e. 

We claim that both of the conditions ^ imply that 7 intersects £. First, if a; € £ 
then this is immediate. Second, if x ^ M(r + e) then we can not have 7([0, £]) C 
M \ £, since this implies that £(7) = 2(7) < r(y) + e, whence a contradiction 
x S M(r + e). Thus 7 intersects £. Let ti > be the smallest i e [0,£] such that 
7(t) € as. Moreover, let f > be the largest t € [0,*i] such that j(t) e M(r). 

We claim that ^ (tI [o,* ] ) > T (y)- First, if to < ^i an d '(7|[o,t l) < T (2/)i then for 
a small t > we have that to + t < t\ and 2(7|[o,t„+t]) < T (y)- Thus 

t + t < t x and j(t +t)eM (r) mt , 

which is a contradiction with the maximality of to- Second, if to = ti and 2(7| ro,t i ) < 
r(y), then y(ti) € M{r) ini n £, which yields a contradiction with Thus 

i(7l[o,to]) > T (2/)- 

We define z := jfa). Then KtI [*o,*i] ) > <M(t),z). Hence 

d(x, 2) + d(z, M(r)) < r(7| [tlifl ) + l(7l[*o,ti]) 

= ^(7l[to,fl) = K7)-K7l[o,to]) 

□ 

Lemma 2. Let VT C M be a neighborhood of M(t) n £. TTien t/iere is ew > 
smc/i t/iat for all e £ (0, ew) 

{x € M (r + e); x satisfies M } C W. 

Proof. We may choose a neighborhood V of M(r) Pl£ such that V C W. Moreover, 
we may choose such ew > that 

(11) e w <d(V,M\W), 

(12) e w <d(M(T,r),j:\V). 

Let e € (0, ew) and let x € M(r + e) satisfy (JsJ) . By Lemma [l] there is z & £ such 
that d(x,z) < ew and d(z, M(r)) < evy. Then zeK, since otherwise (12) would 
be violated. Thus x € W, since otherwise (11) would be violated. □ 

Let ip e C°°(M) be a non-negative function and define the measures 

m^(E) := I ipdm, m^,(E) := / ipdfh. 
J e J E 

The following theorem can be considered as a local version of Theorem [l] 

Theorem 3. Let xo G M(r) D £, and suppose that there is a neighborhood U(xq) 
of xo such that M(r) D U (xq) is an embedded smooth manifold with boundary. 
Then there is a neighborhood V(xq) C M of xq and e(xo) > such that for all 
e e (0,e(x )), 

nv(V(a!o) n M(r + e)) - m^(V(x ) n M(r + e)) 
>mi(e,^,x ) -G(e 2 ), 
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M(r) 




< 5 

Vg) 
Is 



FIGURE 2. A neighborhood of xq and the sets Is and Vs. 



where mi(e,ijj,xo) > 0. Moreover, if iJj(x ) > i/ien 

liminf e _3 ^ 2 mi(e, t/;, xo) > 0. 

Before proving Theorem [3j let us introduce some notation and prove a couple 
of lemmas. Let xq and U(x ) satisfy the assumptions of Theorem [3j and let us 
consider such semi-geodesic coordinates in a neighborhood U C U (xq) of Xq that 
x Q = 0, Sn U = {(x 1 ^ 2 ) € B; x 1 > 0} and 

(13) g(x\x 2 ) = (dar 1 ) 2 + ft(ar 1 ,aj 2 )(da: 2 ) a , 

where B is a neigborhood of the origin in R 2 and h is a stricly positive smooth 
function on B. We may choose Co > such that for all x — (a; 1 , a; 2 ) € J7 and 

y = (y\y 2 ) e 17, 

(14) id(a:,l/) < - z/| + |x 2 - y 2 | < Codfoj/). 



Co 



Let us define 



p{x) :=d{x,M{T)) : ^(x 2 ) :=p(0,x 2 ). 



As M{t)C\U is a smooth manifold with boundary, the map p is smooth in U\M(t) 
near xq, whence also <j) is smooth in a neighborhood of the origin. Non-negativity 
of <p and 0(0) = yield that <j>'(0) = 0. In particular, there is a constant a > such 
that 



(15) 



2\2 



(x 2 ) < a(x 2 ) 



in a neighborhood of origin. Also p(0) = and using boundary normal coordinates 
of U \ Mir) we see that gradp(O) ^ 0. By ^ we have that d x ip(0) > 0. Thus we 
may replace U with a smaller neighborhood of Xq still denoted by U such that p is 
smooth in U, 



(16) 



d x ip(x)>0, x€U\M(t), 



and (15) holds when (0,a; 2 ) G J7. 



Let S € (0, 1). As </>'(0) = 0, there is a neighborhood C K of origin such that 
10' (a? 2 ) I < <VC for all a; 2 € Zj, where C is the constant in (14). We may choose a 
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neighborhood Vg of xo and eg > such that for all x = (x\x 2 ) G M 
B(Vg,eg)cU, B(V s ,e s )ndEc{0}xl s , 

where B(Vg, eg) := {x G M; d(x, Vg) < eg}. Next we study how the set M(r+e)nVg 
stretches in the a^-direction compared to M(r + e)D Vg. We show that the stretch 
is "of magnitude c" in E (Lemma and that it is negligible in M\ E (Lemma [5]). 



Lemma 3. Let e g (0, 6,5) and a; = (x 1 , x 2 ) £VgC\ M(t + e). If x satisfies then 
there is s G Tg such that 

d((x\x 2 ), (0, s))<e- 0(a) < " <K* 2 )) + <^ 2 )- 

In particular, if <fi(x 2 ) > e then | re 1 1 = 0(e 2 ). 

Proof. By Lemma [l] there is z G 9E such that d(x, z) < e — d(z, M(r)). As x € Vg 
and e < £5, we have that z € U. Let us define s G K by (0, s) = z. Then s G 2,5, 
and 

d(x, z)<e- d(z, M(t)) = e - 0(s). 

By (14) we have that \x 2 - s| < C <i(x, z) < C (e - 0(s)) < C e, and by Taylor's 
theorem 

0(x 2 ) = cj>(s) + 0'(s)(x 2 - s) + R{x 2 ,s), 
where |i?(x 2 ,s)| < \\<t>" \\ L ^ C^e 2 = 0(e 2 ). Hence 

e - 0( s ) = e - <j){x 2 ) + (j}'{s){x 2 - s) + R(x 2 , s) 

< e - </>(x 2 ) + -^C (e - 0(a)) + i?(x 2 , s). 



Thus 



( S )< I ^( £ -^ 2 ) + ii: (^ s )) 
= ^(e-#r 2 )) + 0(e 2 ). 



(x 2 ) is 



By Q we have that la; 1 1 < C d((x\ x 2 ), (0, s)) whence l^r 1 1 = 0(e 2 ) if e - 
negative. □ 

Lemma 4. Let e G (0, e«) and x = (a; 1 , a; 2 ) G Vg n M(r + e). J/x G E £/ien 
xl <#^) (e _^ 2))+0(e2) _ 

Proof. If x G <9E then x = 0. Thus we may assume that x G E mt . Let s G 2a be 
as in Lemma [3] Let j(t) — (7 1 (i), 7 2 (i)), i G [0,^], be a path from (0, s) to x in 
the manifold (M,g) that is shorter than e — </>(s). Let a G [0,£) be the parameter 
satisfying j((a,£]) C E m * and 7(a) G <9E. We may assume that j\[ at £\ is a geodesic 
with respect to metric g. 

Note that d(j(t), (0,x 2 )) = 0(e). Indeed, by |l4| we have 

d((0,x 2 ), (0, s)) < C \x 2 -s\< Codx 1 ! + |x 2 - s\) 

<C$d(x,(0,s))<C 2 T('y)=O(e), 
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whence 

d( 7 (t), (0, x 2 )) < d( 7 (t), (0, a)) + d((G, a), (0, x 2 )) 
<T(7) + 0(e)=0(e). 

Moreover, 

c(7(*)) < c(0,x 2 ) + ||gradc|| LOO(Sn[/) d( 7 (f), (0,x 2 )) 
= c(O,x 2 )+0(e), t£(aj], 



whence by ( 13 1 



< 



= (c(0,x 2 ) + O(£)) / |S s 7(s)|^<(c(0,^) + O(e))Z(7) 

< (c(0, x 2 ) + 0(e))(e - 4>(8)) = c(0, x 2 )(e - 0( S )) + G(e 2 ). 
By Lemma [3] we have 

^<c^) (e _^ 2))+0(e2) _ 



□ 



Lemma 5. Let e £ (0, e<s) and a; = (x-^x 2 ) € Vj PI M(r + e). Suppose that 
x <£ M(t + e) U S. TTiera | 1 1 = C(e 2 ). 

Proof. As i ^ M(r + e), we have that 

< d(x, M(t + e)) < p(x) - e. 

Indeed, there is z £ M(T,t), and a path 7 : [0,^] — ► M from 2 to x such that 
J (7) = p{x). Moreover, there is y £ dM such that d(y,z) < r(y). We have 
2(7) > e, since otherwise d(y, x) < d(y, z) + 2(7) < r(y) + e, whence a contradiction 
x G M(t + e). Let < £ [0,1] be the parameter satisfying Z(7|ro,tl) = £■ Then 
d(y,l[t)) < r(y) + e and 7(t) G M(t + e). Moreover, 

d(x,M{r + e))<d(x, 7 (t))<l( 7 \ M ) = l( 1 )-l( 1 \ m ) = p(x)-e. 



By (16) the map, s i-> p(s,x 2 ), is increasing. Note that x 1 < as x ^ E, whence 
e < Kx 1 ^ 2 ) < p(0,x 2 ) = 0(x 2 ). 
Thus Ix 1 ! = 0(e 2 ) by Lemma [3| □ 
Proof of Theorem^ Let e G (0, eg). We define 

V(e) := y 5 n M(t + e) H E, 7(e) := V5 n M(t + e) n E, 
VK(e) := V« n M(r + e) \ E, W(e) := V s n M(r + e) \ E. 
As m^(E) = m^,{E) for a measurable set Ed M\E, Lemma [5] yields 
nv(W(e)) - m*(W(e)) = m^(W(e) \ W(e)) = 0{e 2 ). 
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Let us denote a := (1 — S)^ 1 . By Lemma[4] 

(17) V(e) C {(x 1 , x 2 ); < x 1 < ac(0, x 2 )(e - (f>{x 2 )) + 0(e 2 ), a; 2 G 

We define X a (e) ={i 2 elj; e - 0(x 2 ) > 0} and 

V x {e) := V(e) n {(x\x 2 ) G Va; x 2 G 2^(e)}. 

Then m^(V(e)\^i(e)) = 0(e 2 ). 

Note, that ip( x ) '■— VK^Is^)! 1 ^ 2 is smooth, whence 



li 



ip(x) < ip(0,x 2 ) + 



_ Ce, x G We). 

L°°(sn(7) 



Similarly c(x) 2 V>(ai) < c 2 (0, x 2 )'0(O, x 2 ) + 0(e). We denote c = mings c(x) and 



mi(e»:= / V(0, x 2 )(e - 0(x 2 ))dx 2 



Then 



m^(Vi(e)) = / c(x) 2 ip(x)dx 

r ac(0,x 2 )(t-<j>(x 2 ))+O(t 2 ) 



< 



I 5 (e) JO 



c(x) ip{x)dx dx 



< 



fc(0, x 2 )- 2 i/>(0, + 0(e)) «c(0, x 2 )(e - </>(x 2 ))dx 2 + 0(e 2 ) 
< ac _1 mi(e,ip) + 0(e 2 ). 



By ( 13 1 we have 

{(x\ x 2 ); < x 1 < e - 0(x 2 ), x 2 G l S (e)} C V(e), 

whence 



m^(V(e)) = / ^(:r)da; > / / iptydx^-dx 2 

JV(e) JX e (e)Jo 

= m 1 (e,^)-0(e 2 ). 
By combining the above results we have 

m^(7i n M(r + e)) - ^(Vj n M(r + e)) 

= m^(T(e)) - m*(V(e)) + ^(W(e)) - ro*(W(e)) 
= m v ,(y(e))-m^(y 1 (e))-0( e 2 ) 
> (1 -ac- 1 )mi(e,V) -C(e 2 ). 



As c > 1, we may choose small enough <5 G (0, 1) so that 1 



> 0. We define 



V(x ):=V s , e(x ):=e s , mi(e, ip, x ) := (1 - ac l )mi(e, ip). 



By (15) we have 1(e) := [—ye/a, y'e/a] C ^(e) for small enough e > 0. Note 
that 



C 



V>(0,x 2 ) > V(0) - 



r e, x 2 GX(e). 
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Hence 

mi(e,ip)> _^(0,x 2 )(e-a(x 2 ) 2 )dx 2 

> (^(Q)-O(Vi)) (e-a(x 2 ) 2 )dx 2 

= (V7(0) - O(^)) ^ £ 3/2 = 4^0) £ 3/ 2 _ 0(e 2 } 

Hence -0(0) > implies that 

liminf e _3 / 2 mi(e, ?/;) > 0. 

□ 

Proof of Theorem^ As M(r) n E is compact, there are 

Xi, . . . , xn <E M(t) n E 

such that with the neighborhoods Vj := V(xj), j = 1, . . . ,N, given by Theorem[3j 
we have 



N 



M(r)nEc \JV r 

3 = 1 

Let us choose a neig hborhood of M{t) n E such that W C U£=i V?- Let e > 
satisfy 

e < £vf and e < e(ajj), j = 1, . . . ,N, 
where ew is the constant of Lemma [2j Then by Lemma [2] 

M(t + e) \ W = M(r +e)\lf CM\S. 

Indeed, we have always M(r + e) C M(t + e) and if x € M(r + e) \ then by 
Lemma [2] we have a; S A/(r + e) and x G Af \ E. Hence 

m(M(r + e) \ = m(M(r + e) \ W). 

Let (V'j)jLi be a partition of unity of W subordinate to (Vj)^ =1 satisfying 
Tpj{xj) ^ 0, j — 1, ...,N. Then by Theorem [3] there are ci > and e± > 
such that for all e G (0, ex) 

m(M(r + e) n W) - m(M(r + e) n W) 

N 

= E ( TO ^ (M(t + e) n ty) ~ fh^ (M(r + e) n VJ)) 

> ^rrnfe^.ij) -0(e 2 ) > c l£ 3 / 2 - 0(e 2 ). 
i=i 

Thus for small e > 

m(M(r + e))~m(M(r + e)) _ 1/2 

£ 3/2 - Cl 

□ 
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Figure 3. The domain of influence M(t) is the grey area, that 
is, the upper half plane. The inclusion £ is the union of the lower 
half plane and the triangle of height h. Let g be the Euclidean 
metric and let c(x) = 2, x E E. Then M(r) contains the half disk 
of radius h. 

Example 1. Let us consider the geometry described in Figure^ We denote by T 
the triangle in the figure and suppose that its height h = 1. Moreover, let us denote 
by H the half disk in the figure. Then for small a > we have that m(M(r)) > 
m(M(r)). 

Indeed, we have for small a > that 
1 



i(M(r)nS) 



4 m(M(r) HE) > -(m(T)+m(U)) 
= - (tana + tt/2) > tana = m(T) = m(M(r) n E). 



Moreover, it is always true that 

m(M(r) \ E) = m(M(r) \E) > m(M(r) \ E). 

In the example, E is non-smooth but it is clear that we may smoothen E so that 
the change in the volume is negligible. Thus the example shows that if we allow 
M(t) to penetrate deep into E then the volume m(M(r)) may become larger than 
the volume m(M(r)). 



3. Distance to the inclusion 
To simplify the notation let us denote for r 6 C(dM), 

V(t) := m(M(r)), V(t) := m(M(r)). 
In this section we prove the following two theorems. 

Theorem 4. Suppose that (M,g) is a complete manifold without boundary. Let 
x £ dM, r > and define T r (y) = r — d(x,y). Then 

d(x, E) = max{r > 0; V(r r ) = V(r r )}. 
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Theorem 5. Suppose that (M,g) is simple. Let x £ dM, < h < r and define 
T r,h(x) '■= r and T Tt h(y) '•= h for y =/= x. Then d(x, E) is the supremum of the 
numbers r > with the following property: 

(C) There is h £ (0, r) so that df. V(r r j,) exists for almost all h in [0, ho]. 

By an argument similar with the argument after Theorem [TJ we can choose 
such r £ C(dM) that M{t Ti h) = M(r). Thus the volume data can be used 
to determine if property (C) holds. The claims (i) and (ii) in Theorem [2] are 
consequences of Theorems [5] and |4j respectively. 



Lemma 6. Let (M,g), x £ dM and r r be as in Theorem^ Then B(x,r) D M — 
M{r r ). Moreover, if 

(18) M(r r ) nS^0 and M(r r ) n E m * = 0, 

then M(r r ) is an embedded smooth manifold with boundary near E. 

Proof. We have for all y, z € M that d(y, z) < z) since M C M. Let z € M(r r ). 
Then there is y S 9M such that d(y, z) < r r (y), whence d(x, z) < d(x, y) + d(y, z) < 
r and z £ B(x, r). 

Conversely, let z £ B(x, r) n M. If z £ dM then 

d(z, z) — < r — d(x, z) = r r (z), 

whence z £ M(r r ). Thus we may suppose that z £ M mt . Let 7 : [0, £] — > M be a 
shortest path from x to z. Then there is s £ [0,£) such that 7(s) £ dM and 7|( s ^i 
is a geodesic of (M mt ,g). Let us denote the length of 7 by £(7). Then 

d(7(s), z) < l{-i [sA ) = l(j[ s ,t]) = - Kl\[o,s]) 
<r- d(x,j(s)) = T r (j(s)). 

Hence z £ M(r r ). We have shown that B(x,r) = M(r r ). 

Let us now assume that ( 18 1 holds. Then d(x, S) = r and a point Zq in B(x, r)nE 
is a nearest point to x on the boundary 9E. Hence there is a neighborhood U Zn C M 
of z such that for all z £ U Zo 

z h- > d(x, z) is smooth and V z d(x, z) ^ 0, 

see e.g. the proof of [27} Lem. 2.15]. As dB(x, r) is a level set of z n- d(x, z) we have 
that it is a smooth embedded manifold near zq. This is true for all zo £ B{x, r) n E, 
whence dB(x, r) is a smooth embedded manifold near E. □ 

Proof of Theorem^ If r < d(x, E) then Lemma [(3] yields that 

M(r r ) = B(x, r) n M C M \ E, 

whence also M(r r ) — M (r r ) and m = m on M(r r ). 

If r = d(x,E) then (18) is satisfied and Lemma ^ and Theorem [l] yield that for 
small enough e > 0, V(x r + £ ) < V(-r r+e ). □ 

Lemma 7. Let r ri h be as in Theorem [5| Suppose that (M, g) is simple and let 
ro > 0. Then there is ho > such that r i-> V(r r is smooth near ro /or almost 
all h £ (0, /i ). 
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Proof. We denote 

B(r) := {z G M ; d(z,x) < r], M(h) := {z G M; d(z,dM) < h}, 

where x G <9M is as in Theorem [5j Moreover, we denote by v(y) G T y (M) the 
interior unit normal vector of dM at y and define 

d+S x (M) := {£ € T X (M); \£\ g = 1, (£,v(a:)) fl > 0}, 

tv(0 := max{i > 0; exp x (t£) G M}, £ G d+S x (M). 

As (M, g) is simple we have that tm is smooth and the map 

(19) (t,0 -> exp^tf). < t < r M (0> £ € 9+^(M), 

gives coordinates in M \ {a:}. Moreover, for y = (7m (£)>£) m \ i 2 -} we nave 
(«(») ) ftexp 1B (rii f (00)»<0. Thus 

d t d{ey^> x {tC),dM)\ t=TM{i) = (g^ad y d(y,dM),d t exp x (T M (00)g 

= {v(y),dtexp x (T M (€)£))g + o. 

We define S := {£ G d+S^M); tm(£) > r o/2}- By the implicit function theorem 
and compactness of S there are ho > 0, a finite number of points £j G S, j = 1, . . . N, 
neighborhoods 

and smooth functions dj such that for all £ G Uj, t G Ij and /i G [0, ho) 

d(exp x (t£), dM) = h if and only if t — (Jj(£, h), 

and that the sets Ij x Uj cover the points (tm (£)'£); £ G S. Furthermore, we may 
require that /io < ro/2 and that the boundary normal coordinates are well-defined 
on M(h ). 

Let us define for h G (0, fto) and £ G d + S x (M), 

F(h,£) := exp y (hv(y))\ y={TM{i)t0 . 

Then F gives coordinates in M(ho) \ B(ro/2) and the range of dF at z in dB(ro) H 
M(ho) mt is the whole tangent space T 2 (M). Thus the transversality theorem, see 
e.g. 20, Thm. 3.2.7], yields that F h (£) := F(h,£) is transversal to 8B{r ) n 
M(h ) int for almost all h G (0, ho). 

Let /i G (0, ho) be such that Fh is transversal to dB(ro) fl M(ho) mt . To sim- 
plify the notation we denote <r j (£) := <7,-(£,/j) and e^(£) := d^aj(^). Using the 
normal coordinates (19) we may identify T( t ^(M) = spanjdt, g\}. The sub- 
space T( t ^(dB(ro)) is then spanned by d%. By transversality the vector dF^d^ 
has nonzero dt component at the points (ro,0 G dB(ro) fl M(ho) mt satisfying 
d((ro,£),dM) = h. Moreover, for (t,£) G Uj x Ij the <9t component is 

dt(dF h d £ ) = a'jit). 

Hence o-'j(£) ^ for all £ G 5 satisfying o-j(£) = r . 

Let (4>j(£))jLi, be a partition of unity in S subordinate to the cover {Uj)^ =l . 
The solutions of (7j(£) = r$, £ G supp(0j), can not have an accumulation point 
£oo since otherwise we would have a contradiction Oj(£oo) = 0. Thus the number 
of solutions is finite, and we may choose a finer partition of unity (ipj)fL 1 such 
that for any j — 1, . . . , M there is at most one solution <Tj(£) = r , £ G supp(ipj). 
Moreover, we choose (V'jOjli so that <7j(£) = r implies £ G supp(i/'j) mt . Here we 
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have reindexed the functions aj, allowing aj — Uk for j ^ k, so that aj is defined 
in the support of ipj . 

Let r > r /2. If (t,£) € M \ B(r) then £ G S. Moreover, 

(t,0€M\(B(r)UM(h)) 

if and only if there is j = 1, . . . , M such that £ € supp(^j) and r < t < crj(£). We 
define Sj(r) := {£ € supp(V>j); r < Then 

m(S(r) U M{h)) = m(M) - m(M \ (B(r) U M(/i))) 
(20) =m(M)-J2 m \g(s,0\ 1/2 d S d^. 

J = 1 JS 3 (r) Jr 

As <7y(£) does not vanish in the solution set crj(£) = r we get easily that Sj(r ) = 
implies Sj(r) = for r near ro- Moreover, Sj(ro) — supp(V'j) implies S'j(r) = 
supp(^ij-) for r near ro. We see that the terms in the sum (20) corresponding these 
two cases are smooth. 

Let us choose coordinates (—1, l)9a4(6 d + S x (M) and let us consider such 
j = 1, . . . , M that there is a solution a?o € supp(^j) mt of aj (a) = ro. By the inverse 
function theorem aj 1 is a smooth function near ro. Let us denote p := aj 1 and 
define for small e > 0, 

U := {(r, a) G (r - e, r + e) x (-1, 1); p(r) < a}. 

The set U is open. Let us assume for a moment that p'(ro) > 0. Then (r,a) € {/ 
implies that r < <jj (a) and the function 

(r,a) i-> ^-(a) / |g(s, a)| 1/2 ds 



is smooth on i7. 

Let us show for a smooth function 6 that 

V(r) := / 6(r, a) da 

JS 3 (r) 

is smooth near r . We have for small t > and r near r that 

Sy(r) \ Sj(r + ()={ae supp(^); r < (Jj(a) < r + 1} 
= {p{r),p{r + t% 

whence 

V(r) -V(r + t) = / b(r,a)da- / b(r + t,a)da 

= j b(r,a)da- / b{r,a)da-t / 9 r &(r, a)dQ; + C(t 2 ) 

JSj(r) JSj(r+t) JSj(r+t) 

rp(r+t) 

= b(r, p(r))(p(r + t) - p(r)) + / 0(t)da 

Jp(r) 

-t l d r b(r,a)da + 0{t 2 ). 

Js 3 {r+t) 
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The case t < is analogous, whence for r near ro, 

V*(r) = -b(r, p{r))p'{r) + [ d r b(r, a)da. 

Js 3 (r) 

Clearly, the first term is a smooth function of r near r . The second term is 
differentiable by the above argument. Hence V is smooth near ro by induction. 

The case p'(ro) < is analogous and we see that r i->- m(B(r)UM(h)) is smooth 
near ro. □ 

Proof of Theorem^ Let r > satisfy r < d(x, E). Then for h < <i(E,9M) and 
r near ro, we have V^tv^) = V(r r /j). Hence condition (C) holds by Lemma 7 

Let r = d(x, E) and < /i< d(E, 0M). Then (|6]) is satisfied for L = 9M and 
T — T ra .h- Moreover, using the coordinates (19) we see that M(t) is an embedded 
smooth manifold with boundary near E. By Lemma[7]there is < ho < d(E,dM) 
such that r H ► V(r rj ^) is smooth near ro for almost all ft, € [0, ho]. Let us fix such 
h G [0,h ] and denote V(r) := V{r r . h ) and V"(r) := V(r r , h ). Then 

y(r + e) - 2T/(r ) + V>(r - e) 
e 2 

^ V^o + e) - V>o + e) V(r + e) - 2V/(r ) + V(r - e) 
e 2 e 2 
The second term converges to V"(ro) and the first term diverges by Theorem [I] 
Hence V"(ro) does not exist, and we see that (C) does not hold. □ 

4. The direct problem 

In this section we study the regularity of the solution u* of the wave equation 
Q. In particular, we see by combining results [35] with results [M], that the 
Neumann-to-Dirichlet operator A2T is well defined on L 2 ((0,2T) x dM). 

We define a quadratic form by 

Qg(u 7 v) :— / (du,dv)gdm, u, v G H (M), 
Jm 

where du denotes the exterior derivative of u and (•, •)§ is the inner product on the 
contangent bundle given by ~g. Then 

Ag : H\M) -> (H^M))* 

with homogeneous Neumann boundary conditions is defined by 

(AgW,-y) L 2 (M;ffl) = -Qg(u,v), u,v € if^M). 

By [Ml Thm 9.5] the solution of the problem 

(21) dfu-A g u = FeL 2 (0,2T-(H 1 - e (M)Y), u t<Q = 0, 

satisfies u G C([0, 2T]; H e {M)), where G (0, 1) and u t<0 = stands for u| t=0 = 
d t u\ t= o = 0. Let 6 G (0, 1/2), / G L 2 ((0,T) x 9M) and define 



(22) (F(t),v):= f(t,x)v\ dM (x)dS g , u e ^"'(M), 

./9M 

where dS'g is the Riemannian volume measure of (ffM, ffl ajy )' Then _F is in the 
space L 2 (0, 2T; (H 1 ~ e (M))*) and the solution u — u f of (21) solves the equation 
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QSJ) in variational sense, see equation (9.21a) in [34], p. 288]. In particular, the map 
/ i y w (T) is a compact operator 

L 2 ((0,T) x 5M) -> i 2 (M). 

The Laplace-Beltrami operator A ff of (M, 5) has smooth coefficients, and by [32 
the solution of the Neumann problem 

d 2 v - A g v = 0, v t<0 = 0, d v v = / G L 2 ((0, 2T) x 9M) 

satishcs zj 6 C([0, 2T]; iJ Q (M)) for any a < 3/5. Let us choose a cutoff function 
X G C°°(M) such that \ — 1 m a neighborhood of 9A/ and that g = 5 in the 
support of x- Let us define w := w — Then w satisfies (21) with 

F = (8 2 - A-sX-xv) = - X {d 2 t - A> - [A ff , x ]v = [ X , A>. 

Note that the commutator [ X , A 9 ] is a differential operator of order one. More- 
over, as X — 1 i n a neighborhood of 9M we have that [x, A ff ] = in the same 
neighborhood. Thus there are x € Cq 3 (M) and C > such that 

([x,A») £/(M)xCoo(M) = ([x,A>,x>)^, (M)xC ao (M) 

= [X, A g ]*x4>) < M\h"{m) \\\x,&g\*x4>\\ H -°>(M) 

< c \\v\\ Ha ^ M ) ||^||hi-«(m) • 

Hence F= [x,A 9 ]wg C([0, 2T]; (H l ~ a (M))*) and to G C([0, 2T]; H a (M)). Choos- 
ing a € (1/2, 3/5) we see that A 2 t is a compact operator 

L 2 ((0,T) x DM) -> i 2 ((0,T) x 0M). 

Note that, if / G Q?°((0, 2T) x M) then u and F are smooth and by 34, Thm 
8.2] 

(23) u f €C([0,2T];H\M)), d t u f G C([0, 2T]; L 2 (M)), 
t V, Aju' G L 2 (0,2T; (iJ x (M))*). 

5. Computation of the volumes 

In this section we show that the minimization algorithm of |37j works also with 
a piecewise smooth metric tensor. We denote 

(24) Rf(t):=f(2T-t), Jf(t):=^J^l L (t,s)f(s)ds, 

K := JA 2T - RA 2T RJ, If(t) := 1 (0 ,t)(*) / /(*)<*», 

Jo 

where L := {(*, s) e R 2 -,t + s < 2T, s > i > 0}. Moreover, we denote by 1+ the 
adjoint of / in 



S : = L 2 ((0, 2T) x dM; dt <g> dS, 



9)1 



where dS g is the Riemannian volume measure of the manifold (dM, g\dM)- 

Let us denote by (•,•) and ||-|| the inner product and the norm of S. For r G 
C(dM) we define S(t) to be the set of / G S satisfying 

(25) supp(/) c {(t,y) G [0,T] x dM; t G [T - r(y),T]}. 



INVERSE OBSTACLE PROBLEM FOR THE WAVE EQUATION 



19 



We claim that the regularized minimization problem, 

(26) axgmin ((/,#/) - 2(//, 1) + a ||/|| 2 ) , a > 0, 

has unique solution f a , a > 0, that is also the solution of the linear equation, 

(27) (P T KP T + a)f = P T I + l, a > 0, 

where P T is the orthogonal projection from S to S(t). The proof is based on the 
following two reformulations of the Blagovestchenskii identity, 

(28) (u'(T), u h {T)) L 2 (M . dfh) = (/, Kh), 

(29) {u f {T),l) L , {M , dK) = (//,l), 

where and it' 1 are the solutions of ^ corresponding to the boundary sources 
f,h e L 2 ((0, 2T) x DM), respectively. For the original formulation of the identity 
see [5] . The equations ( 28 ) and ( 29 1 together with compactness of / i— > w (T) imply 
the existence and uniqueness for the minimization problem ( 26 ) by using exactly 
the same argument as in |37j . 

Let us verify next the equations ( 28 ) and ( 29 1 . By continuity of the map / i->- 



v,f(T) and of the operators A 2 t, J, R and / it is enough to show (28) and ( |29| ) for f 
and h in Cq°((0, 2T) x dM). Then and u h have the regularity properties (23). 
We have 

{d 2 t -d 2 s ){ui{t),u h {s)) LWfn) 

= (5 t 2 M / (t),u' l (s)) (fl - 1(M)) » xffl(M) - (" / ( i ):^ U ' l ( S ))/ i -l(M)x(ffl(M))* 

f{t,x)u h \ dM {s 1 x)dS g 



+ q(u h (s),uf(t)) 



<)M 



h(s,x)u f \ dM (t,x)dS g 



(/(*> x)A 2T h(s, x) - h(s, x)A 2T f{t, x)) dS g , 



whence we obtain ( 28 ) by solving a one dimensional wave equation with vanishing 
initial and boundary conditions. Note that, the right hand side of this wave equa- 
tion is in L 2 ((0, 2T) x (0, 2T)), whence the solution is in C([0, 2T];H 1 (0, 2T)) and 
w(T,T) is well defined. This follows also from ([23]). Moreover, 



dUu f (t),l) L 2 (M . dfh) = (dju^t),!) 



(H 1 {M))*xH 1 (M) 



-q(uf(t),l)+ f(t,x)dS : 



DM 



f(t,x)dS g 



and we obtain ( 29 1 by solving an ordinary differential equation with vanishing initial 
conditions. 



We claim that for all r E C T (dM), 



(30) 



lim u* a (T) = 1 j 



in L 2 (M), 



where f a , a > 0, are the solutions of (|27| ). See Figure [4] for examples of compu- 
tational such solutions of (|27|) that u^ a (T) approximates ljjjV T )- Note that A 2 t 
determines the volume data(|5j) by ( |30| and (28), whence we have proved Theo- 
rem [2] after proving (30). Moreover, the convergence (30) can be proved exactly in 
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Figure 4. Waves w (T) approximately one on M(r). The bound- 
ary source / has been obtained by solving a discretized version of 



(27) with the conjugate gradient method. Here g is the Euclidean 
metric tensor and we have chosen t G C(dM) so that M(t) is a 
half plane. On left, there is no inclusion. In middle, S is a disk 
and c(x) = 5 for x G S. On right, E is the same disk as in the 
middle but c(x) = 1/5 for x G S. 

the same way as the corresponding result in [ 37] given the following approximate 
controllability result. 

Lemma 8. Lei T C dM be open, T > and de/ine r(y) := Tl r (y), y G <9M. Then 
£/ie embedding 

{u f (T); f G C °°((0,T) x T)} C L 2 (M(r)) 

is dense. 

This result is well-known when the metric tensor is smooth, see e.g. [57]. We 
outline the proof in the appendix below. 

Appendix A. Unique continuation 

In this appendix we employ the unique continuation result by Tataru |47j to 
prove Lemma [8] We begin by proving two other lemmas. 

Lemma 9. Let T > 0, U C M mt be open and define V := U D <9£. Let 

u G C([—T, T]; H 1 ^)) H C 1 ([— T, T];L 2 {U)) 

be a solution of d^u — A-gu = 0. Then the traces u \^_t T)xr an< ^ ® uU \~(-TT)xr 
of n|(-T,T)xc/ns vanish if and only if the traces u\^_ T T ) xr and 9i,u\7_ TT \ xT of 
u\(-T,T)xU\X vanish. 

Proof. Let us assume that the traces u\^_t T)xr an< ^ ^ u l(-TT)xr an d vanish 
on (— T, T) x r. We denote by the distance function of the Riemannian manifold 
(T,,g). Let S G (0, T). By Tataru's unique continuation, see [UJ Thm. 3.16] for a 
formulation suitable for our needs, we have that u(t) vanish on 

V s :={xeUnT,; d E (x,T) < 6} 

for all t G (— T + 5,T — 5). Thus it satisfies the wave equation 

(31) d;u-A g u = Q on {-T + S,T- 5) x V s U (Z7\ s )> 
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Again by unique continuation there is a neighborhood Us of T in M such that u 
vanish in (— T + 26, T — 25) x Us for small 5 > 0. Hence the traces u \7_x t)xT and 
d u u\7_ T T w r vanish on (— T, T) x T, and we have shown the necessity part in the 
claim of the lemma. The sufficiency can be shown analogously after choosing a 
smooth Riemannian metric on M that is an extension of g. □ 

Lemma 10. Let T > and let 

u e Cd-T^^H^M)) n C 1 ([-T,T];L 2 (M)) 
be a solution of d 2 u — Ajit = satisfying for an open T C dM 

w|(-T,T)xT = 9 l/ u|(_ T ,T)xr = 0. 

Then u(t, x) vanish when 

(t, x) E (— T, T) x M and d(x, T) + \t\ < T. 

Proof. Let (t ,x ) e (— T, T) x M satisfy d(x ,T) <T — \t \. Then there is y G T 
and a shortest path 7 : [0, •£]—>• M from ?/ to Xo such that ^(7) < T — \to\. The open 
set 7 _1 (Af \ S) \ {0,^} is a union of open disjoint intervals (aj,bj), j € J, where 
either J = N or J = {1, . . . , N} for some N e N. 

Let us show that J = N leads to a contradiction. As the intervals are disjoint, 
we have 

00 

Thus Z(7|[a j ,6 j ]) -> as j -t 00. Moreover, j(aj),j(bj) e <9£, and 7l(aj,6 3 ) is a 
geodesic in M \ E. Let us consider the semi-geodesic coordinates (13) in a neigh- 
borhood U of 7(aj). For large enough j we have that j([aj,bj}) C U. Let us define 
7j(t) := (0,7 2 (t)), £ € [aj,6j]. Note that jj is a smooth path from 7(07-) to j(bj). 

We will show that JCtIcc^)) > J(7j) f° r large j, which is a contradiction since 
7 is a shortest path. By ( 13 ) we have 



J a j 



Moreover, for £ G [oj, 6j] 

7 1 (£) - d( 7 (t),0E) < d( 7 (*),7(ai)) < K7l[a J)t ]) <T(7|[, J)6j 
Thus 7 1 roo , n — ^ as j — > 00, whence for large j 

II II L (a,j,bj) 

VH0,l 2 (t)) < min c(z) y/H-yit)), t g [a, ,b 3 ]. 



zeds 



Moreover, £ 3 \da 2 {t)\dt > since otherwise 7 2 (£) is constant for £ G [aj,bj] and 
7| [a •,&•] is a loop. Hence for large j we have the contradiction 



c(0,7 2 (£)) 



We have shown that J is finite. 
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We may renumber the intervals (dj, bj), j = 1, . . . , N, so that 

ai < bi < a 2 < b 2 < ■ ■ ■ < cin < 6jv- 

Then we may repeatedly use unique continuation together with Lemma [9] either in 
(M \ or in (T,,g) to see that u(t) vanish near 7(5) for \t\ < T ~ ^ (tI [o,s] ) - I n 
particular, u(to) vanish near j(£) = xq since |to| < T — □ 



Given Lemma 10 the proof of Lemma[8|is almost identical with the proof of [271 
Thm. 3.10]. Note that the boundary of M (Tip) is of measure zero by the proof of 
(371 Lem. 21. 
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